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Abstract 

. In this paper, we calculate the branching ratios and CP-violating asymmetries for — > 

^ I {p^ , oj , (t))r]^'^ decays in the perturbative QCD (pQCD) factorization approach. Numerically we 

I found that (a) the pQCD predictions for the CP-averaged branching ratios are Br{Bs p^-q) ~ 

■ 0.07 X 10^*^ , Br{Bs p^r]') « 0.10 x 10"^ , Br{Bs coi]) ^ 0.02 x 10"^ Br{Bs cori') ^ 

^ '■ 0.13 X 10-^ Br{Bs (t^rf) « 2.7 x 10"^ and Br{Bs H) ~ 2.0 x 10"^; (b) the gluonic 

contributions are small in size: less than 3% for B — > {p, uj, (f))rj decays, and about 10% for 
B {p,uj,(j))ri' decays; and (c) the pQCD predictions for the CP-violating asymmetries of the 
considered decays are generally not large in magnitude. The above predictions can be tested in 
^ ' the forthcoming LHC-b experiments at CERN. 
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I. INTRODUCTION 



The theoretical calculations and experimental measurements of the rare B meson de- 
cays play an important role in testing the standard model (SM), probing CP violation 
of B meson system and searching for possible new physics beyond the SM. At present, 
about 1000 million events of B meson pair productions and decays have been collected by 
BaBar and Belle collaborations. In the forthcoming LHC experiments, a huge amount of 
B meson events, say around 10^^ ~ 10^^, are expected, and therefore the rare B meson 
decays with a branching ratio around 10~^ can be observed with good precision. Another 
advantage is that the heavier Bg and Be mesons and b-baryons, besides the Bu and Bd, 
can also be produced and studied at LHC [H. 

By employing the generalized factorization approach^, 0] or the QCD factorization 
(QCDF) approach [3], about forty Bg — > /ii/i2 {hi stand for light pseudo-scalar or vector 
mesons ) decay modes have been studied in the framework of SM jl, 0, 0, [1] or in some 
new physics models beyond the SM [qI]. In this paper, we will study the Bg — > p'^r]^'^ , 
ujT]^''^ and (pT]^'"^ decays in the pQCD factorization approach. In principle, the physics for 
the Bg two-body hadronic decays is very similar to that for the B^, meson except that the 
spectator d quark is replaced by the s quark. 

For Bg — > Co", 0)?7*^'^ decays, the Bg meson is heavy, setting at rest and decaying 
into two light mesons (i.e. p° and r]^'^ ) with large momenta. Therefore the light final 
state mesons are moving very fast in the rest frame of Bg meson. In this case, the short 
distance hard process dominates the decay amplitude. We assume that the soft final state 
interaction is not important for such decays. The smallness of FSI effects for B meson 
decays into two light final state mesons has been put forward by Bjorken based 
on the color transparency argument [llj], and also supported by further renormalization 
group analysis of soft gluon exchanges among initial and final state mesons [l2|. With the 
Sudakov resummation, we can include the leading double logarithms for all loop diagrams, 
in association with the soft contribution. Unlike the usual factorization approach, the 
hard part of the pQCD approach consists of six quarks rather than four. We thus call it 
six-quark operators or six-quark effective theory. Applying the six-quark effective theory 
to Bg meson decays, we need meson wave functions for the hadronization of quarks into 
mesons. All the collinear dynamics are included in the meson wave functions. 

This paper is organized as follows. In Sec. [Tll we give a brief review for the PQCD 
factorization approach. In Sec. Illlt we calculate analytically the related Feynman dia- 
grams and present the decay amplitudes for the studied decay modes. In Sec. HV] we 
show the numerical results for the branching ratios and CP asymmetries ofBg — > p'^r]^'^ 
, UT]^'^ and (pT]^'^ decays and comparing them with the results obtained in the other two 
methods mentioned above. The summary and some discussions are included in the final 
section. 



II. THEORETICAL FRAMEWORK 



The pQCD factorization approach has been developed and applied in the non-leptonic 
B meson decays for some time [ll|, [l^, 14, 15]. This app roach is based on fc^ factor- 
In this approach, the decay 



This a 

ization scheme, where three energy scales are involved [1 
amplitude is factorized into the convolution of the meson's light-cone wave functions, the 
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hard scattering kernels and the Wilson coefficients, which stand for the soft ($), hard(H), 
and harder(C) dynamics respectively. The hard dynamics (H) describes the four quark 
operator and the spectator quark connected by a hard gluon. This hard part is charac- 
terized by a/AMb^, and can be calculated perturbatively in pQCD approach. The harder 
dynamics (C)is from mw scale to m^^ scale described by renormalization group equation 
for the four quark operators. The dynamics below a/AMb^ is soft, which is described by 
the meson wave functions ($). While the function (H) depends on the processes consid- 
ered, the wave function is independent of the specific processes. Using the wave functions 
determined from other well measured processes, one can make quantitative predictions 
here. Based on this factorization, the decay amplitude can be written as the following 

A{Bs ^ M1M2) ~ Jd%d%d% TT[CmBAki)^MAk2)^M,{h)H{kuk2,h,t)], (1) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C{t) is Wilson coefficient of the four quark operator which 
results from the radiative corrections at short distance. The functions $m and H are 
meson wave functions and the hard part respectively. 

For the Bg meson decays, since the b quark is rather heavy we consider the Bg meson 
at rest for simplicity. It is convenient to use light-cone coordinate {p~^,p~ , Pt) to describe 
the meson's momenta, 

= ±P^), and Pt = ip\p^). (2) 

Using these coordinates the Bg meson and the two final state meson momenta can be 
written as: 

Mr Mr n Mr 

Pb. = ^{1,1,0t), P, = -^{l,r%Or), = -^(0, 1 - rj, Ot), (3) 

respectively, where Vp = nip/MB^ and the light pseudoscalar meson masses have been 
neglected (here we just take the decay Bg — >• p^t] as an example). Putting the light (anti-) 
quark momenta in Bg, and 7] mesons as fci, ^2, and ^3, respectively, we can choose 

ki = (xiP+, 0, kiT), k2 = {X2P^, 0, k2T), k3 = (0, 0:3^3", k3T). (4) 

Then the integration over k^ , k^ , and k^ in eq.([T]) will lead to 



A{Bg ^ p rf) ^ J dxidx2dxzhidhih2dh2hzdh?j 

■Tr [CmBX^iM)^p{x2M)%{xzM)H{xiM,t)St{x,)e-'^'^] , (5) 

where hi is the conjugate space coordinate of kiT-, and t is the largest energy scale in 
function H{xi, bi,t), as a function in terms of Xi and 6j. The large logarithms {In miy/t) 
coming from QCD radiative corrections to four quark operators are included in the Wilson 
coefficients C{t). The large double logarithms (In^Xj) on the longitudinal direction are 



summed by the threshold resummation 15| , and they lead to St{xi) which smears the end- 
point singularities on Xi. The last term, e~^^^\ is the Sudakov form factor resulting from 
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overlap of soft and coUinear divergences, which suppresses the soft dynamics effectively. 
Thus it makes the perturbative calculation of the hard part H applicable at intermediate 
scale, i.e., Mb^ scale. We will calculate analytically the function H{xi,bi,t) for our six 
decays in the first order in as expansion and give the convoluted amplitudes in next 
section. 

For the considered Bg — > p^r]^'\ujr]^'^ and (pr]^'^ decays, the low energy weak effective 



Hamiltonian H^jf for the b ^ q transition with q = {d, s) can be written as [16 

r r 10 ■ 

We specify below the operators in Tie// for 6 — > transition: 

Oi = daYLup ■ u^jf^Lba , O2 = daj^Lua ■ u^j^Lbp , 

O3 = d^YLba ■ Eg' Q'lB^f^Lq'p , = d^l^Lbp ■ Y.q' Qp^f^Wa > 

O5 = d^YLba ■ Eg' Q-'pli^RO-'d ' = daYLbp ■ Eg' Q-'pl i^R(l'c > 

O7 = IdaYLb^ ■ Eg' eq'q'p^f^Rq'f^ , Og = ^d^YLbp ■ Eg' eg' q'f^-f f,Rq'^ 

O9 = IdaYLb^ ■ Eg' dq'q'p'^^Lq'p , do = ^d^^Lbp ■ Y.q' dq'Q'plt^Lq'^ 



(6) 



(7) 



a 1 
I 

la 1 



where a and (3 are the SU{?>) color indices; L,R = (1 =F75) are the left- and right-handed 
projection operators. The sum over q' runs over the quark fields that are active at the 
scale fi = 0{mf,). Since we here work at the leading twist approximation and leading 
double logarithm summation, we will also use the leading order (LO) expressions of the 
Wilson coefficients Ci(u) (i = 1, . . . , 10), although the next-to-leading order Ci{fi) already 
exist in the literature [l6| . This is the consistent way to cancel the explicit n dependence 
in the theoretical formulae. 

For the renormalization group evolution of the Wilson coefficients from higher scale to 
lower scale, we use the formulae as given in Ref. [3] directly. At the high mw scale, the 



leading order Wilson coefficients Cj(Mvv) are simple and can be found easily in Ref. [16 
In pQCD approach, the scale 't' may be larger or smaller than the rrib scale. For the case 
of nib < t < mw, we evaluate the Wilson coefficients at t scale using leading logarithm 
running equations, as given in Eq.(Cl) of Ref. [ij]. For the case of t < rrib, we then 
evaluate the Wilson coefficients at t scale by using Ci{mh) as input and the formulae 
given in Appendix D of Ref. [3] • 

For the wave function of the heavy Bs meson, we take 

= + MBj750B,(ki). (8) 

Here only the contribution of Lorentz structure (pBsO^i) is taken into account, since the 
contribution of the second Lorentz structure (pB^ is numerically small and can be neglected. 
The distribution amplitude (pBs iii Eq. ([8]) will be given lately in Eq. (147|) . 

For B — >■ Vt]^'^ decays, the vector meson V = (p, 0, to) is longitudinally polarized. The 
relevant lon gitu dinal polarized component of the wave function for p meson, for example, 
is given as [17|, 

*P = y= y [mpM^) + ]^pKi^)] + mp(Pp{x) } , (9) 
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where the first term is the leading twist wave function (twist-2), while the second and 
third terms are sub-leading twist (twist-3) wave functions. For the case of V = u and 
their wave functions are the same in structure as that defined in Eq. (Q, but with different 
distribution amplitudes. One can find the distribution amplitudes 0^,,^ and (j\f^{x) in next 
section. 

For 1] and r]' meson, the wave function for dd components of r]^'^ meson are given as 



18 



^75 

V2K 



Vdd 





- V ■ n)(p^jx) 



(10) 



where P and x are the momentum and the momentum fraction of rj^ respectively, while 
and represent the axial vector, pseudoscalar and tensor components of the 



'7dd' ^'7dd 



wave function respectively, and will be given explicitly in next section. Following Ref. [18 



we here also assume that the wave function of rj^^ is the same as tt wave function based on 
SU(3) flavor symmetry. The parameter C, is either +1 or —1 depending on the assignment 
of the momentum fraction x. 

Before we proceed to do the perturbative calculations, we firstly give a brief discussion 
about the cj) — uj mixing , as well as the rj — rj' mixing and the gluonic component of the 
7]^'^ mesons. 

For the vector (p — u meson system, we choose the ideal mixing scheme between 0(1020) 
and u{782) 



LU 



—j= [uu + dd) 
V 2 



— ss. 



-uu 



+ dd)/^/2 [19|. 



since the current data support this ideal mixing scheme p^. The quark contents of p 
meson is chosen as p° = 



For the 77 — 77' systern, there exist two popular mixing basis: the octet-singlet basis and 
the quark- flavor basis [13, 21]. Here we use the quark-flavor basis lO] and deflne 

riq = {uu + dd)/ V2, r]s = ss. (12) 
The physical states t] and t]' are related to rig and ris through a single mixing angle 0, 



f/(0) 



Vs 



cos ( 
sin ( 



— sm I 

cos 



Vg 

Vs 



The corresponding decay constants fg, fs, f^''^ and /^/** have been deflned in Ref. [20 



(13) 



as 



< 0\sYl5s\r]s{P) > 

< 0\qYl5q\v'-'\P) > 

< 1 S7^75S I r/W(P) > 



~^frf') P^ ' 



(14) 



(15) 
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while the decay constants Z^'** and f^l are related to fq and fs via the same mixing matrix, 



fq fs 

J Tj J rj 

fl fs 

Jr)' Jr}' 



u{4>) 



fs 



(16) 



The three input parameters fq, and in the quark-flavor basis have been extracted 



from various related experiments [20|, l21 



/, = (1.07 ± 0.02)/^, /, = (1.34 ± 0.06)A, 



39.3° ±1.0°, 



(17) 



where /^r = 130 MeV. In the numerical calculations, we will use these mixing parameters 
as inputs. 

As shown in Eq. ( |T3l) . the physical states rj and rj' are generally considered as a linear 
combination of light quark pairs mm, dd and ss. But it should be noted that the rj' 
meson may have a gluonic component. Following Ref. 221], we also estimated the possible 
gluonic contributions to B ^ (p, uj, (f))r]^'^ decays induced by the gluonic corrections to 
the B — > 1]^'^ transition form factors 2^ and found that these corrections to both the 



branching ratios and CP violating asymmetries are indeed small: less than 10%. 

Frankly speaking, on the other hand, we currently still do not know how to calculate 
reliably the gluonic contributions to the B meson decays involving t]' meson as final state 
particle. For the studied decay modes in this paper, we firstly consider only the dominant 
contributions from the quark contents of r]^'^ meson, and then take the subdominant 
contribution from the possible gluonic content of r]^'^ meson as a source of theoretical 
uncertainties. 



III. PERTURBATIVE CALCULATIONS 

In this section, we will calculate the hard par t Hjt) for the considered decays. Following 
the same procedure as being used in Refs. |23l . |24J . |25| . |26 | and taking B — > p'^rj^'^ decays 



as an example, we will calculate and show the analytical results for all decay amplitudes, 
and then extend the results oi B p^r]^'^ decays to other decay modes under study. 

As illustrated in Fig. [H there are eight type Feynman diagrams contributing to the 
Bg p^T]^'^ decays. We first calculate the usual factorizable diagrams (a) and (b). Opera- 
tors Oi, O2, O3, O4, O9, and Oio are {V — A){V — A) currents, the sum of their amplitudes 
can be written as: 

rl poo 

F^r, = -SttCfM^J'p / dxidxs / hidhih^dh2,(t)Bs{xiM) 

Jo Jo 

• { [(1 + X3)0^(x3, h) + r^(l - 2x3)(0j'(x3, 63) + ^^(a^s, h))] 
as{tl)heixi, X3, 61, 63) exp[-5'ab(tg)] 

+2r^(p^{x3,bs)as{tl)heix3,Xi,bs,bi)exp[-Sab{tl)]} . (18) 

where C^? = 4/3 is a color factor and = m'^' /Mb^ or = itIq'^'^ /Mb^ since depends 
on the quark components in t]. The function h^, the scales and the Sudakov factors Sab 
are displayed in the appendix. From diagrams Fig. 1(a) and 1(b), one can also extract 
out the form factor Fq''^^"". 
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(h) 



FIG. 1: Diagrams contributing to the Bs — > p^i] decay (diagram (a) and (b) contribute to the 
Bs rjss form factor F^'^'''~). 

The operators O5, Oq, O7, and Og have a structure of (V — A)(y + A). In some decay 
channels, some of these operators contribute to the decay amplitude in a factorizable way. 
Since only the vector part of (ViA) current contribute to the scaler meson production, 

{p\V + A\B){r]\V - A\0) = {p\V - A\B){r]\V - A\0), (19) 

the result of these operators is the same as Eq. f[T5]l . and therefore we find easily that 

FS = Ferj . (20) 

In some other cases, one needs to do Fierz transformation for these operators to get right 
flavor and color structure for factorization to work. In this case, we get (5" + P){S — P) 
operators from {V — A){y + A) ones. For these {S + P){S — P) operators. Fig. 1(a) and 
1(b) will give 

F^^ = . (21) 

For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 63 can be performed by using 5 function 5(63 — 61), leaving 



7 



only integration of 61 and 62- For the {V — A){V — A) operators, the corresponding decay 
amphtude is: 



32 

Mf^n = —^TvCpMfj / dxidx2dx3 / bidbib2db2 
V6 Jo Jo 



■ <j)Bs{xi, bi)(t)p{x2, h2){x'i [(t)^{xz, bi) - 2rr,(j)^{x3, bi)] 

■ Ois{tf)hf{xi, X2, X3, bi, 62) exp[-Scd{tf)] } ■ (22) 

For the {S — A)(V + A) operators, the decay amphtudes read 

64 f°° 
= —TrCpM^rp / dxidx2dx3 / 6id6i62C?&2 ^^^(a^i, &i) 
V6 "Jo Jo 

■ { [X2(I)^{X3, bi) {(j)pix2, &2) - 0p(.«2, ^2)) + {ix2 + X3) (^^(Xs, 61) 

• <Ppix2, 62) + 0^(0:3, 61)0* (a;2, 62)) + (3:3 - X2){(j)^{x3, 61)0* (2:2, 62) 
+0j(x3, bi)(j)p{x2, 62)))] as{tf)hf{xi, X2, X3, bi, 62) exp[-5'cd(t/)] }, (23) 
= -Me, . (24) 

For the non-factorizable annihilation diagrams 1(e) and 1(f), again all three wave func- 
tions are involved. Here we have all three kinds of contributions. M^^ is the contribution 
containing operator type {V — A){V — A), while M^^ and M^j^ is the contribution con- 
taining operator type {V — A){V + A) and {S — P){S + P), respectively 

32 f°° 

Mar^ = -='kCfM% I dxidx2dx3 / bidbib2db2(j)Bs{xi,bi) 

v6 'Jo Jo 

■ {- [X2(t)pix2, b2)(f)^{x3, 62) + Tprr, {{X2 - X3) {(p^{x3, 62)0p(a;2, ^2) 
+(f>'^{x3, b2)(pp{x2, 62)) + (2 + 0:2 + X3)(l>^ix3, 62)0^(2^2, ^2) + ("2 + 

X2 + X3)0^(x3, 62)0^(2^2, ^2))] as{tj)h'j{xi, X2, X3, bi, 62) exp[-5'e/(t/)] 

+ [•^■30p(a;2, b2)4>^{x3, 62) + TpTr^ {{X3 - X2) {(j)^ {xs, b2)(t)\{x2, 62) + 
0J(X3, b2)(f)l{x2, &2)) + (2^2 + X3) {(j)^{x3, b2)(l)p{x2, 62) + 0^(2^3, ^2) 

• (I>l(x2, 62)))] q;,(4)/i}(xi, X2, X3, 61, 62) exp[->Se/(4)] } . (25) 



32 /* /'CO 

M^i = -—ttCfM^^ J dxidx2dx3 J bidbib2db2(t)BX^iM) 

■ { [rp{2 - X2)(f)^{x3, 62) {(f>l{x2, b2) + 0p(x2, 62)) + r,(a;3 - 2)0^(^2, 62) 

• (0^(2^3, h) + 0^(2^3, ^2))] Q;s(iJ)/iJ(xi, X2, X3, 61, 62) exp[-5'e/(ij)] 
+ [x2rp(t)^{x3, &2) (05,(2:2, &2) + (t>^p{x2, ^2)) - a:3r,0p(2;2, ^2) 

• (0j' (X3, 62) + <(X3, ^2))] a.(t})/i}(xi, X2, X3, 61, &2) exp[-5e/(t})] } . (26) 
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32 f°° 

M^^ = -^ttCfM^^ / dxidx2dx3 / hidhih2dh2 4>Bs{^i-,^i) 

V 6 Jo Jo 

• { [x30p(x2, 62)0^ (a;3, ^2) + rprn(l)^{xs, 62) ((a;2 + 2:3 + 2)0^(x2, 62) 
-(x2 - 0:3)0* (0:2, &2)) + rprr,(f)'^{xs, 62) ((a;3 - X2)(l)pix2, 62) + (a;2 

+X3 - 2)0^(a;2, 62))] X2, 0:3, 61, 62) exp[-S'e/(tj)] 

+ [i-X2)(l)pix2, 62)0^ (a^s, ^2) + ^^p?^,, {{xs - X2) (05'(a;3, &2)0p(a;2, ^2) 
+0j'(a;3, &2)0p(a;2, ^2)) - (a;2 + 2:3) (0j'(a;3, &2)0p(a;2, ^2) + 0j'(a;3, ^2) 

0p(a;2,&2)))] as(t})/i}(xi,X2,X3,6i,62)exp[-S'e/(t})]} . (27) 

The factorizable annihilation diagrams 1(g) and 1(h) involve only p and rj wave func- 
tions. Again decay amplitude Fa-^ is for {V — A)iy — A) type operators, while the F^^ 
and F^^ come from the {V — A)(y + A) and {S — P){S + P) type operators, respectively 

I* 1 /"OO 

Far^ = MCpfBs^h, / '^a;2rfx3 / b2db2hdbs 
Jo Jo 

■ { [X30p(x2, h2)(t)^{x^, h) + 2rpr^0^(x2, 62) ((1 + X3)0J'(X3, 63) 

+ (x3 - 1)0^(^3, 63))] asitl)ha{x2, X3, 62, ^s) exp [-S",,?, (t^)] 
- [a;20p(a;2, 62)0^(2:3, ^a) + 2rpr^(f)^ {x^, 63) ((1 + a;2)0p(a;2, 62) 

+ (x2-l)0p(x2,62))]as(te)/ia(x3,X2,63,&2)exp[-53/,«)]} , (28) 
-^aTT^ ~ ~Farii (29) 

/•l /"OO 

f5 = -IQ-nCpfEMk / ^^2^X3 / b2db2hdh 

Jo Jo 

■ { [2rp0p(a;2, &2)0;^(a;3, ^a) + rr,Xs(f)p{x2, ^2) (0j'(a;a, ^a) 

-0j'(a;a, &a))] as(te)^a(a;2, 2:3, &2, &a) exp[-Sghitl)] 

+ [x2rp(l)^ix3, 63) (0p(a;2, 62) - 0p(a;2, ^2)) + 2r^0p(x2, 62) 

• 0^(^3,63)] as(tg)/ia(a;3,a;2,&3,&2)exp[-S'gh(te)]} • (30) 

When we exchange the position of p and r]^'^ mesons in Fig. [H then only 4 annihila- 
tion diagrams l(e)-l(h) can contribute to -B — > p^r]^'^ decays. The corresponding decay 
amplitudes from the non- factorizable annihilation diagrams 1(e) and 1(f), where it is the 
p° meson who picks up the spectator s quark, can be written as 

32 f°° 

Map = -=ttCfM^ / dxidx2dx'i / 6ic/6i62'^&2 0_Bs(^1) ^1) 

V 6 'Jo Jo 

■ { [X30p(x3, b2)(t)^{x2, 62) + TpT^ ((^3 - X2) {(j)^ ix2, &2)0p(a;3, h) 

-F0j'(a;2, &2)0p(a;3, ^2)) + (3^2 + 3:3) (0j'(a;2, &2)0p(2;a, h) + 0j'(a;2, ^2) 

■ 0p(a;a,&2)))] as(t})/i}(a;i, 0:2, X3, 61, 62) exp[-5'e/(t})] - [a;20p(a;3, 62) 
■0;^(a;2, ^2) + ((^2 - x-i) (0^(a;2, &2)0p(a;a, ^2) + 0j'(a;2, ^2) 

■ 0p(x3, 62))) + rpTr, ((2 + X2 + X3)0j'(x2, &2)0p(a;3, h) - (2 - X2 - X3) 

■ 0^(2^2, &2)0p(a;3, ^2))] as(t'j)h'j{xi, X2, X3, 61, 62) exp[-S'e/(tj)] } , (31) 
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M^p = —^ttCfM^^ / dxidx2dx3 / hidbib2dh2 4>Bs{^ii^i) 
v6 ^ Jo io 

■ { [a;2r^0p(x3, 62) (0j'(a;2, h) + 0j'(a;2, ^2)) - x^^Vp {(f)p{xs, 62) 
+0p(a;3, ^2)) <P^ix2, 62)] X2, 0:3, 61, 62) exp[-5'e/(t})] 

+ [(2 - X2)r^0p(x3, 62) «(a;2, ^2) + <(a;2, &2)) - (2 - a:3)rp (0^(0:3, 62) 
+0* (x3, 62)) 0;^(a;2, &2)] X2, X3, 61, 62) exp[-Sef(tj)] } , (32) 



ap 



32 



dxidx2dx3 / 6i(i6i62'i&2 0_Bs(^i? ^1) 

■ { [a;20p(x3, h2)(t)^{x2, ^2) + ^pr^ ((^2 + X3) (0^(x2, b2)(j)pixs, 62) 

+0J(X2, 62)0p(a;3, &2)) + ix2 - X3) (0J'(X2, 62)05,(2:3, h) + (t>^{x2, 62) 

■ 0^(x3,62)))] as(t})/i}(xi,X2,X3,6i,62)exp[-S'e/(t})] - [a;30p(x3, 62) 
• 0;^(a;2, 62) + ^^p^r, ((a;3 - X2) (0j'(x2, 62)0p(a;3, 62) + 0j'(a;2, 62)0p(a;3, 62))) 
+rpr^ ((2 + X2 + X3)0^(x2, 62)0p(a;3, 62) + {X2 + x^ - 2)0^(X2, 62) 

■ 0p(a;3, 62))] a5(t/)/t/(a;i, X2, 0:3, 61, 62) exp[-5'e/(t^)] } . (33) 

For the factorizable annihilation diagrams 1(g) and 1(h) after the exchange of and 



7]^'' mesons, the corresponding decay amphtudes can be obtained directly through the 
links with their counterparts F^^, and F^^ 



ap 



-F 

aril ap 



aril 



pP2 ^ pP2 
ap ar) ' 



(34) 



Combining the contributions from different diagrams, the total decay amplitude for 
Bg p^rj decay can be written as: 



er) 



e. ( Ci + ic2 



eri 



iuC2 — 6 I ^2^^ 



C7- 



ap y 



^mC'2 — it 2^10 



i^2(0) 

Fi(0)+(m5+M5) 



i^2(0) 



-it 



+ (i^a^ + F, 



ap J 



iu[C, + -C2 



3 13 1 



Fi( 



(35) 



where ^„ = V^J^Kd, 6 = V'^feV'td, and Fi(0) = cos 0/2 and ^2(0) = - sin 0/^2 are the 
mixing factors. The Wilson coefficients Ci should be calculated at the appropriate scale 
t by using the formulas as given in the Appendices of Ref. 271 . 
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Following the same steps, one can derive the total decay amplitude for Bg oorj decay: 



er) 



1 



2^ 1^ 1^ 1^ 
3 2 ^ 6 2 ^ 6 



i^2(0) 



C2 ~ C,t { 2C4 — 2Cq — -Cs + 2^10 



C2 — \ 2C4 + 2^10 



i^l(0) 



+ (m5+m2) 



-6 ( 2C6 + -Q 



2 21111 

— 'Ct 1 2C3 + -C4 — 2C5 — -Cq — -C7 — -Cs + -Cg + gC'io 



Fi{4>). (36) 



The individual decay amplitudes [F^rj, M^^, ■ ■ ■ ) in Eq. fl36l) can be obtained easily from 
those as given for the case of Bg p^rj decay by the simple replacements 



/p 



ful^ f 



(37) 



Note that the difference in the quark components of the two vector mesons p° and uj has 
been taken into account. 

For Bs 0?7*-'-* decays, all eight Feynman diagrams when the rj^'^ or meson picks up 
the spectator s quark. For the later case, we firstly consider the factorizable diagrams 1(a) 
and 1(b). The decay amplitude induced by inserting the {V — A){y + A) operators is 

/* 1 noo 

Fecp = —SnCpM^^ / dxidx^ / bidbib^dbs (f)B^{xi,bi) 
" Jo Jo 

■ { [(1 + X3)M^3, bs) + r^(l - 2x3){(f)l{x3, 63) + (P'^xs, 63))] 
asit].)he{xi, X3, bi, 63) exp[-Sab(t].)] 

+2r^(f)''^{x3,b3)as{tl)he{x3,xi,b3,bi)exp[-Sab(tl)]} , (38) 

where = m^/MB^- One can also extract out the Bg ^ (j) form factor A^"'^'^ from the 
diagrams 1(a) and 1(b) when cj) meson picks up the spectator s quark. For other diagrams, 
the relevant decay amplitudes can be written as 



(39) 



/»1 /»oo 

— IQ-nCpM'jjjri / dxidx^ / 61C/6163C/63 0b^(xi, 61) 
Jo Jo 

■ { [00(a;3, ^3) + ((X3 + 2)0^(X3, 63) - X30*^(X3, 63))] 

■as{tl)K{xi, x^, 61, 63) exp[-Sab{tl)] + (xicp^ixs, 63) 

+2r<^0^(x3, 63)) a^(tg)/le(x3, Xi, 63, 61) exp[-S'afe(te)]} 



(40) 
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32 f°° 

Mecf, = —j=T:CfMq^ I dXidX2dX3 / 61(^6162^62 

v6 "Jo Jo 

&l)0^(a;2, 62) {X3 [00(X3, 61) - 2r^0*^(x3, 61)] 

■ a,(t/)/i/(xi,X2, 0:3,61, 62) exp[-S'cd(t/)]} , 
= 0, m5 = M,^. 

Finally, the total decay amplitude of Bg (pi] decay can be written as 



(41) 
(42) 



e« ( Ci + -C2 



'^t { 2C3 + -C4 — 2C5 — -Cq — -C7 — -Cg + -Cg + gC'io 



/si^2(0) 



4 4 1112 2 

^'3 3 ^3 2^6 3 3 



e« (C2) - 6 ( 2C4 + 2^6 + + iCio 



-6 ( C3 + C4 + a - ^Cg - - iCio^ i^2(0) 

"■Ct ( C3 + C4 — Cq + -Cs — -Cg — 2^10 



-6 ( C, - -C, 



+ V2Merj 

+v^(ms + m5) 
+v^(m5 + m5) 

+y2(F,^ + F„ 



i^2(0) 



i^2(0) 



"6 { C3 + C4 — -Cg — 2^10 



i^2((/') 



4 4 

— C3 -|- — C*^ 
3 3 



+ ^C, + \c, - - ^Cio 



i^2(0) 



"6 ( 2^*5 + Cq — -C7 — -Cs 



i^2(0). 



i^2(0) 



(43) 



The individual decay amplitudes {F^r^, Mg^, ■ ■ ■ ) in Eq. (H3i) can be obtained easily from 
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those as given for the case of Bg — > p^t] decay by the simple replacements 

fp — ^ /<^, /J — ^ /J, rup — > 
0, 0^, 0* 0*^, 0^ 0; (44) 

The total decay amplitudes for Bs P^v\ ^^v' 4>v' can be obtained 

easily from Eqs. (!35l) . (!36|) and (H3|) by the following simple replacements: 

Fi(0) = COS0 — ^ F{(0) = sin0, 
F2(0) = -sin0/V2 — > F!^{<p) = cos 0/^2. (45) 

Note that the possible gluonic component of r]' meson has been neglected here. We will 
estimate its effects in next section. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



A. Input parameters and wave functions 

We use the following input parameters in the numerical calculations 



A 



(^4) 
MS 





mn 



250MeV, fp 

1.4GeV, fB, = 
0.782GeV, 
1.02GeV, f4, = 
0.770GeV, M. 



205MeV, 

: 230MeV, 
= 195MeV, 
: 237MeV, 



leOMeV, 



f 

= 130MeV, 
/J = 140MeV, 



220MeV, 



5.37GeV, M, 



w 



80.42GeV. 



(46) 



For the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, here we adopt the 
Wolfenstein parametrization for the CKM matrix up to 0{X^) with the parameters 
A = 0.2272, A = 0.818, p = 0.221 and f] = 0.34o[l9|. 

Bg meson is different from B meson due to the heavier strange quark (compare to u,d 
quark) which induces the SU(3) symmetry-breaking effect. This effect is considered to be 
small and the distribution amplitude of Bs meson (given in the following formula) should 
be similar to that of the B meson, 



5)5^ (x, 6) = Nb,x^{1 - xfexp 



Ml 



2uo 



B. 



(47) 



where ub^ is a free parameter in nature. After considering the constraints from some Bg 
non-leptonic decays [2^, we here use ub^ = 0.55 ± 0.05 in the numerical calculation. 

For the light meson wave function, we neglect the h dependant part, which is not 
important in numerical analysis. We choose the wave function of p{uj) meson similar to 
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the pion case 2^: 



= [l + 0.18C|(2x-l) 

0pm(^) = ^ {3(2a; - 1)^ + 0.3(2a; - 1)^ [5(2x - 1)^ - 3] 



+0.21 [3 - 30(2x - if + 35(2x - 1)^]} , 

- 2x) [I + 0.76(10a;2 - 10a; + 1)] . 



For the wave function of meson, we use (29|: 



J 6 



2\/6 
4^6 



3(1 - 2x)2 + 1.68Q (1 - 2x) + 0.69 



1 + (1 - 2x) In 



1 — X 



3(1 - 2x)(4.5 - 11.2x + 11.2x2) + 1.38 In 



X 



1 — X 



where the Gegenbauer polynomials are defined by 

c'ko = li^e-i), 

= ^ (35^^ - 30e' + 3) . 



(48) 

(49) 
(50) 

(51) 
(52) 
(53) 

(54) 
(55) 



For 1] meson's wave function, and represent the axial vector, pseudoscalar 

and tensor components of the wave function respectively, for which we utilize the results 
from the light-cone sum rule 30| including twist-3 contribution: 



- ^) {l + [5(1 - 2xr - 1] 

+a^"^ [21(1 - 2x)^ - 14(1 - 2xf + l] | , 

^-IpL) [3(1-2x)^-1] 
- ^pL«^") [35(1 - 2x)^ - 30(1 - 2xr + 3] | 



2v^ 



/, 1 + - 30n. 



+ - ( -3r]3UJ3 
-^=U1 - 2x) 



20 ^''dd 



7 



6 



20 '^''dd s^'Jdd 



(56) 



with the updated Gegenbauer moments [31 



0.115, al^^ = -0.015, 
ttIt, / ttIq'^'^ , = 0.015 , 



-3.0 



(57) 
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We assume that the wave function of uu is same as the wave function of dd. For the 
wave function of the ss components, we also use the same form as dd but with m^'"' and 

fs instead of mg'^'^ and /g, respectively. 



B. Branching ratios 



For the decays we have considered here, the decay amplitudes in Eqs. ( !35l) , fl36|l and 
( |43|) can be rewritten as 

M = V:,VusT-V:,VtsP = V:,VusT[l + ze'^^^^'^] 

where 

P 



(5^ 



T 



(59) 

is the ratio of penguin to tree contributions, 7 = arg — Yr'^ll is the weak phase (one of 

the three CKM angles), and 5 is the relative strong phase between tree (T) and penguin 
(P) diagrams. 

From Eq. fl58|) . it is easy to write the decay amplitude for the corresponding charge 
conjugated decay mode 

M = Vuiy:sT - Vt,V:,P = Vu,V:,T [l + ze^^-^^'^] . (60) 

Therefore the CP-averaged branching ratio is 

Br = {\M\^ + \M\^)/2 = \VubV:jf [1 + 2^ cos 7 cos (5 + ^2] , (61) 

where the ratio z and the strong phase 6 have been defined in Eqs.( l58l) and ( l59l) . 

Using the wave functions and the input parameters as specified in previous sections, it 
is straightforward to calculate the CP averaged branching ratios for the considered decays: 

(62) 

(63) 



Br{B, - 


- P%) = 


[0.07: 


h0.03/, , 
-0.02 l^Bs 


<+0.05 
'-0.03 






X 


10 


-6 


Br{Bs - 


> p'i) = 


[0.10: 


h0.04/ , 


<+0.07 
'-0.04 


[rris 




X 


10 


-6 


Br{B, - 


-^urf) = 


[0.21: 


h0.04/ 
-0.03 l^Bs 


v+0.14 
'+0.00 


[rris 


) ±0.00( 


7)] 


X 


10 


Br{B, - 




[0.13: 


h0.04/ 


,+0.05/ 




X 


10 


-6 


Br{Bs - 


= 


[2.66: 


hi. 10/ , 


v+1.45 
'-0.78 


[rris 




X 


10 


-5 


Br{B, - 




[2.00: 


h0.80/, , 


,+1.42 
'-0.73 






X 


10 


-5 

7 



(64) 



where the main errors are induced by the uncertainty of uje^ 
130 ± 30 MeV and 7 = 60° ± 20°, respectively. 



Br{Bs 
Br{B, 

Br{B, 
Br{B, 

Br{B, 
BriB^ 



0.55 ± 0.05 GeV, m. 



1 Ref. 0: 










- p'v) = 


[0.122- 


0.151] 


X 


10 


- p'v') = 


[0.123- 


0.160] 


X 


10 


s ^UJT]) = 


[0.006 - 


0.025] 


X 


10 


^UT]') = 


[0.010 - 


0.075] 


X 


10 


s^(j)T]) = 


[0.088 - 


0.417] 


X 


10 


^ 0^0 = 


[0.024 - 


0.149] 


X 


10 



(65) 
(66) 
(67) 
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y (degree) 



7 (degree) 



FIG. 2: The CP averaged branching ratio (in unit of 10 ^) of Bs p^rj decay as a function of 
CKM angle 7. (a) is for = 130 MeV and lob, = 0.50 GeV (dotted curve), 0.55 GeV(sohd 
curve) and 0.60 GeV(dashed curve); and (b) for ujb^ = 0.55 GeV, and rus = 100 Mev(dotted 
curve), 130 MeV (sohd curve) and 160 MeV (dashed curve). 




- -I 0.04 r " 

30 60 90 120 150 ISO 30 60 90 120 150 ISO 

y (degree) y (degree) 

FIG. 3: The same as Fig. [2] but for Bg — > p^r]' decay. 



One can see that the pQCD predictions are basically consistent with the corresponding 
QCDF predictions for the first four decays, but much larger than the QCDF predictions 
for Bs —>■ <f)i] and (pr]' decays. The difference is about two orders! Further study about 
this difference is under way. 

In Figs. we present the PQCD predictions of the CP averaged branching ratios of 
the considered decays and their dependence on the variations of ub^, rris and the CKM 
angle 7. 
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FIG. 4: The same as Fig. [2] but for Bs — > ujij decay. 



m 
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FIG. 5: The same as Fig. [2] but for Bg cor]' decay. 



C. CP- violating asymmetries 

Now we turn to the evaluation of the CP-violating asymmetries of the considered 
decays in PQCD approach. For these neutral decay modes, the effects of — mixing 
should be considered. 

For _B° meson decays, the CP-violating asymmetry of — > fcp decay is time 
dependent and can be defined as 



Acp = 



r (^°(At) fcp) - T (i?°(At) ^ fcp) 
T (4o(At) ^ fcp) + F (i?o(At) ^ fcp) 
Afp cos(Am,At) + A^]? sin(Am,At), 



(68) 



where Am,; is the mass difference between the two i?^ mass eigenstates. At = tcp — ttag 
is the time difference between the tagged B^ (5°) and the accompanying B^ (5°) with 
opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct and 
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20.0 

-e- 



CD 



FIG. 6: The same as Fig. [2] but for Bs (prj decay. 
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FIG. 7: The same as Fig. [2] but for Bg — > cpr]' decay. 



mixing induced CP-violating asymmetries A'^p and A^^p can be written as 



A 



dir 
CP 



lA 



CP\ 



1 



12' 



1 + |Acp. 

where the CP- violating parameter Xcp is 

VtbVts{fcp\Heff\B^) 



A mix 
^CP 



2Im{X 



CP) 



1+ A 



CP I 



12' 



A 



CP 



o2n 



l + ze 



i{S--y) 



Vt,v:Afcp\Heffm 



1 + 2e*(^+T) ' 



(69) 



(70) 



Here the ratio z and the strong phase 5 have been defined previously. In PQCD approach, 
since both z and 5 are calculable, it is easy to find the numerical values of and A^j? 
for the considered decay processes. 

The PQCD predictions for the direct and mixing-induced CP-violating asymmetries 
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are: 



^^p(5. - p'v') = [-25.31^:«(c.bJ^8^:^k)1^:^(7)] x lO"^ (7i) 

- cry) = [7.6l^i(a^i.Jli^iK)^^i(7)] X 10-^ 
A'l^^ABs - a;r/') = [-7.61° j(a;Bjl|^K)l?;3(^)] x lO'l (72) 

^^p(5. - 0V) = [-o.5l°:^(c.BjlS:iK)i°i(7)] X lo-^ (73) 
^cp (5. - p%) = [-4.2l!|(.;,,jli°:?K)iJ:^(7)] X lo-^ 

A— (5. ^ p%') = [18.7+°:^KJli^K)l^i(7)] X 10-^ (74) 

A— (E.^cr^) = [8.0+°:?(c.Bj:i^(m.)l^:}(7)] X 10-^ 

A— (5. ^ ^V) = [24.6lJ;^(^^J+nK)^6:8(7)] X 10"^. (75) 

^CP - 0r^) = [-0.5l°:°(u;^Jl°:l(m.)l°:l(7)] x 10-^ 

A— (5, ^ 0V) = [0.5l°:^(^BjlS:i(m.)l°:}(7)] x 10-^ (76) 

where the dominant errors come from the variation of = 0.55 ± 0.05, rris = 130 ± 30 
Mev, 7 = 60° ± 20°. One can see from these numerical results that the CP- violating 
asymmetries of the considered decay modes are generally not large in size. This is rather 
different with the cases of Bd and Bu meson decays, where the CP-violating asymmetries 
are generally large in magnitude. 

As a simple comparison, we show the QCDF predictions for the direct CP-violating 
asymmetries of the six considered decays as given in Ref. 0]: 

-16.9 X 10^2 





-A) 


A^UBs - 




A'SUBs - 




A'^UBs - 


UJT]') 


Al^UBs - 




Afp{B, - 





-33.0 X 10"^ (77) 
0.76 X 10"^ 

4.4 X 10"^ (78) 
23.2 X 10-^ 

-58.3 X lO'l (79) 

Since the relevant measurements are not available at present, and the theoretical un- 
certainties in both pQCD and QCFD factorization approaches are still very large, it is 
too early to make a meaningful comparison now. One has to wait for the starting of 
experimental measurements and improvements of the theoretical calculations. 

In Fig. [6] we show the 7— dependence of the direct and mixing induced CP- violating 
asymmetries for B^ — »• p*^?7^'^ and o;?]^'-' decays. Since the pQCD predictions of the CP 
asymmetries are sensitive to many input parameters, the lines in Fig. [6] should be broad- 
ened accordingly. 



D. Effects of possible gluonic component of r]' 



Up to now, we have not considered the possible contributions to the branching ratios 



and CP-violating asymmetries of B 
component of rj^'^ 18|, l32l . 



(p°, uj, (f))r]^'^ decays induced by the possible gluonic 
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FIG. 8: The direct and mixing-induced CP asymmetry (in percentage) of B'^ p^rj (solid 
curve), p^T]' (dotted curve), torj (dashed-curve) , and ujr]' (dash-dot-dot curve) decay as a function 
of CKM angle 7. 



By using the formulae as given in Ref. [22], we calculate the gluonic contributions to 
B ^ 7] and B ^ rj' form factors. For B {p,uj,(f))ri decays, the gluonic contributions 
are negligibly small, less than 3%. For B — > {p,u!, (f))ri' decays, the gluonic contributions 
to the branching ratios are only about 10%. The central values of the pQCD predictions 
for B — s> (p, 07,0)77' decays after the inclusion of gluonic contributions are the following 



BriB, ^ p%') 
Br{Bs UJT]') 
Br{B, ^ (1)7]') 



[0.12l°;°^(c.^j: 

[o.i41°;°^(..bJ: 
[2.iitS:r7KJ: 



0.07 
0.04 

-0.06 
0.03 

-1.45 
0.75 



(m, 



.+0.00(^)] 



t(7)] 



-0.00 1 



X 10 

X 10" 
X 10" 



-6 



(80) 
(81) 
(82) 



Again, the gluonic contributions are small in size, which is consistent with previous 
results for other B meson decays involving rj or rj' as one or two final state mesons 23 , 



V. SUMMARY 

In this paper, we calculated the branching ratios and CP-violating asymmetries of 
Bg p^r]^'\ Bg ujri^'\ and Bg — > (prf^'^ decays in the pQCD factorization approach. 

Besides the usual factorizable diagrams, the non-factorizable and annihilation dia- 
grams as shown in Fig. ([T]) are also calculated analytically. Although the non-factorizable 
and annihilation contributions are sub-leading for the branching ratios of the considered 
decays, but they are not negligible. Furthermore these diagrams provide the necessary 
strong phase required for a non-zero CP-violating asymmetry for the considered decays. 

From our calculations and phenomenological analysis, we found the following results: 

• For the branching ratios of the six considered decay modes, the pQCD predictions 
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are: 




(83) 



(84) 



(85) 



Here the various errors as specified previously have been added in quadrature. The 
pQCD predictions for the first four decay channels agree well with those obtained 

by employing the QCDF approach. For last two decays, however, there are large 
differences between the pQCD and QCDF approach, which will be tested in the 
forthcoming LHC-b experiments. 

• The gluonic contributions are small in size: less than 3% for B (p, 4>)Tj decays, 
and about 10% for B — > {p,u, (j))rj' decays. 

• For the CP-violating asymmetries, the pQCD predictions are generally not large in 
size, but still have large theoretical uncertainties. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function /ij's, coming from the Fourier transformations of H^^\ 



ha{x2, X3, 62, ^3) = Kq {i^Jx2XzmBM) " &2)-ft'o (^A/is^^B.&s) -^0 ii^/x^mBM) 



he{xi,X3,bi,b3) 



Ko {^fx^niBM) [Gipl - b3)Ko {^/x^niBM) h iV^mB.b^) 
+9{b3 - bi)Ko {y/x'3mBM) h {\fx^rnBM)] St{x3), 



(Al) 



+61(62 - b3)KQ {i^fx^niBM) h {iV^^B.bs)] St{x3), 



(A2) 




0{b2 - bi)h{MB,y/^bi)Ko{MB^V^b2) 
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h]:{xi,X2,X3,bi,b2) 



+ 



9{bi - b2)Ko{iy/X2X3biMBjloiiy/X2X3b2MB, 

Ko(Mb,F(2)6i), 



for > 0' 

fH«(MB.,^F|^6i), for ^ < 



(A4) 



/lJ(Xi,X2,X3,6i,62) 



6{bi - b2)Ko{iy/x2X3biMB,)loiiy/x2X3b2MB^ 
Ko(Mb,F(3)6i), 



+ (61 ^ 62 



fH«(MB.,/|F/^3J6i), for F'^, < 



for Fj)>0^ 



(A5) 



(3) 



(3) 



where Jq is the Bessel function, Kq and Jq are modified Bessel functions Ko^—ix) 



-{7i/2)Yo{x) +i{7r/2)Jo{x)] h[,^^(z) is the Hankel function, H, 
F(j)'s are defined by 

(Xi - X2)X3 , 
(Xi - X2)X3 , 

Xi + X2 + X3 - X1X3 - X2X3 . 



^(1) 

^(2) 
^(3) 



Jo(^) + z Yo(2;), and 

(A6) 
(A7) 
(A8) 



The threshold resummation form factor St{xi) is adopted from Ref.|17 

2i+2T(3/2 + c 



St{x) 



v^r(l + c) 



[x(l -x)]^ 



(A9) 



where the parameter c = 0.3. This function is normahzed to unity. More information 



about the threshold resummation can be found in reference [13|, [15 
The Sudakov factors used in the text are defined as 



Sab{t) = s (^XimBjV2,bi^ + s (^X3mBjV2,b3^ + s (^{1 - X3)mBjV2,b: 



1 



Scdit) 



-ln(6iA) -ln(63A) 
s (^XimBjV2,bi^ + s {x2mBjV2,b-^ + s [{I 

+s {x3mBjV2, 61 j + s (^(1 - X3)mBjV2, bi 
1 



(AlO) 



2,6. 



Pi 



21nJ^^ + lnJ^^ 
-ln(6iA) + -ln(62A) 



(All) 



Sefit) 



S XiTUb, 



2, 61 j + S {x2'mB, /V2M)+s{{l-X2)mBjV2,h 
+s {x3mBjV2, 62) + s (^(1 - X3)mBjV2, 62) 



/3i 



lnJm + 21n.^^(^/^) 



ln(6iA) 



ln(62A) 



Sgh{t) = s (x2mBjV2,b2^ + s {x3mB^ /V2,b3)+s{{l-X2)mBjV2,h 



+ s{{l-X3)niBjV2M) 



In 



ln(t/A) 
-ln(63A) 



In 



ln(t/A) 
-ln(62A) 



(A12) 



(A13) 
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where the function Sb^i SpO, S^{,) used in the amphtudes are defined as: 



(A14) 



^po(t) = six2P^, h) + s ((1 - X2)P^, h)+2 ^7 {asift)) , (A15) 



(A16) 



where the so called Sudakov factor s{Q,b) resulting from the resummation of double 
logarithms is given as: 



(A17) 



with 



A = Cp— + 



67 TT^ 10 
Y ~ y ~ 27 



nj + -/^oln 



3 7B 



a. 



TX 



2 a,, /e^^^-i 



B = ^ In 

3 TT 



(A18) 
(A19) 



Here 7^; = 0.57722 ■ ■ ■ is the Euler constant, Uf is the active quark fiavor number. For the 
detailed derivation of the Sudakov factors, see Ref.|3^. The hard scale tj's in the above 
equations are chosen as: 

t\ = max(^/x^mij^, 1/61, 1/63) , 



tl = max(v/iTmi3,, 1/61, 1/63) , 
tl = max(v/i^mB^, 1/62, l/fes) , 
= max{^/x^mB,,l/b2,l/h) , 



(A20) 
(A21) 
(A22) 
(A23) 

tf = max(^/x7x^mB^, y/iix^mB^, 1/61, 1/62) , (A24) 
t] = max{^/x2X3mB,, I/&1, 1/&2) , (A25) 
= max(Va;i + 0:2 + X3 - X1X3 - X2XsmB^, y/x2X3mB,, I/&1, 1/62) • (A26) 

They are given as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 
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